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Abstract 

Nonstationary and nonequilibrium processes are considered on the basis of an 
Enskog-Landau kinetic equation using a boundary conditions method. A nonsta- 
tionary solution of this equation is found in the pair collision approximation. This 
solution takes into account explicitly the influence of long-range interactions. New 
terms to the transport coefficients are identified. An application of the boundary 
conditions method to hydrodynamic description of fast processes is discussed. 
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The development of methods to construct a theory for nonequilibrium pro- 
cesses in dense gases, liquids and plasmas is an important direction in the 
modern theoretical physics. Moreover, the construction of kinetic equations 
for such classical and quantum systems still remains to be a major problem 
in the kinetic theory. It is complicated additionally in the case of dense gases, 
liquids and plasmas, where kinetics and hydrodynamics are closely connected 
and should be considered simultaneously [1-5]. 

An approach for construction of kinetic equations from the first principles of 
statistical mechanics, namely from the Liouville equation, has been developed 
in [6,7]. Another approach for obtaining kinetic equations for dense systems, 
which is based on ideas of papers [6,7], has also been proposed [1] and gener- 
alized in [2,3]. Here, the formulation of modified boundary conditions for the 
BBGKY hierarchy is used taking into account corrections connected with lo- 
cal conservation laws. On the basis of this sequential approach, a new Enskog- 
Landau kinetic equation has been obtained for an one-component system of 
charged hard spheres. There is a considerable interest of an application of 
this kinetic equation for description of transport processes in dense systems of 
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charged particles as well as in ion melts and electrolyte solutions. The normal 
solution and transport coefficients for this equation have been found in the pa- 
per [8] using the Chapman-Enskog method. The same approach has been used 
for a many-component system of charged hard spheres in the paper [9] with 
more detailed calculations for a two-component system as well. At the same 
time, as is well known, the Chapman-Enskog method allows one to find the 
transport coefficients in a stationary case only. Similar drawbacks are peculiar 
to the Grad method [10,11] which is oftenly used to solve kinetic equations 
next to the Chapman-Enskog method. 

In this paper, the Enskog-Landau kinetic equation for a system of charged hard 
spheres is investigated. To find the normal solution in a nonstationary case, the 
so-called boundary conditions method is used, which has been introduced in 
[12,13]. As a result, transport coefficients equations in the nonstationary case 
are written. A limiting case of the stationary process is considered. A brief 
comparison of the obtained transport coefficients with those known previously 
from the Chapman-Enskog method is given. 

Let us consider the Enskog-Landau kinetic equation for a one-component sys- 
tem of charged hard spheres [8] : 



where f\{xi;t) is the one-particle distribution function. The right-hand side 
of this equation is the so-called generalized Enskog-Landau collision integral, 
where each term can be considered as a separate collision integral. Their struc- 
ture are as follows: 

Ie ( x u t) is the collision integral of the Enskog theory RET [14]: 



where a is a hard sphere diameter, g denotes a vector of relative velocity 
for two particles, r 12 is a unit vector along the direction between centres of 
particles 1 and 2, 

v' 1 = v 1 + r 12 (r 12 -g), g = v 2 -vi, 
v 2 = v 2 - r 12 (r 12 ■ g), r 12 = |r 12 |~V 12 ; 




^7 + v itt~ \ fi (xii t) = I E (xi, t) + I MF (xi, t) + I L (xi, t) , 



(1) 



I E (xi, t) =a 2 J df 12 dv 2 (r 12 g) {r 12 g) x 
{92 (ri, ri + r l2 a- 1) j\ (n, v[; t) j\ (n + r 12 a, v 2 ; t) - 
92 (n, ri - r 12 a; t) f 1 (r 1 , Vi, t) f 1 (r 1 - r 12 a, v 2 ; t) }, 



(2) 
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lMF(xi',t) is the collision integral of the kinetic mean field theory KMFT 
[15,16]: 



Imf (an; *) = ^ f dx 2 ~^~^"^92 (n, r 2 ; t) f\ (an; t) fi (x 2 ; t) , (3) 

where & 1 (\r 12 \) is a long-range part of the interparticle interaction potential; 
II (xi',t) is generalized Landau collision integral [2,8]: 
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fi(xi,t) fi(x 2 ;t). 



It is necessary to note that the quasiequilibrium binary correlation function g 2 
takes into account the full interaction potential (hard core part plus long-range 
Coulomb tail). 

One of a major problem at the correct derivation and solution of kinetic equa- 
tions is their consistency with local conservation laws of particle density (or 
mass), momentum, total energy and substantiation of hydrodynamic equa- 
tions and incomprehensible calculation of transport coefficients via molecular 
parameters. These conservation laws for classical systems in general have the 
structure as in [17]. 

To find a solution of the Enskog-Landau kinetic equation (1) using one or 
another method, there is necessary to take the advantage of local conserva- 
tion laws in corresponding approximations. So doing the expressions for ki- 
netic coefficients will be defined through densities for momentum flow tensor 
II (r;t) and energy flow vector j e (r;t) on the basis of solution fi(xi,t) and 
corresponding approximations for g 2 (ri, r 2 ; t). As far as we find the solution 
that corresponds to linear hydrodynamical transport processes by gradients of 
thermodynamical parameters, densities of momentum flow tensor n (r;t) and 
energy flow vector jg(r; t) could be determined immediately with the help of 
kinetic equation (1) without general formulas from [17]. To this end it is conve- 
nient similarly to [2], to introduce the following hydrodynamical parameters: 
density n(ri,t) (or mass density p(ri,t)), hydrodynamical velocity V(ri,t) 
and density of kinetic energy u)k{r\\ t). Multiplying initial kinetic equation (1) 
by hydrodynamical parameters and integrating with respect to Vi, one can 
obtain the equations for these parameters in the form: 
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d o(r 1 ;t)=-4-V(r l ;t), (5) 



p(r i; f)df rv 15 ' dr x 
P (n; (r i; f)= :P (n; t) , (6) 

P (ri; t) ^ (n; *)= -^-g (n; t) - P (n; f) : ^-V (n; f) , (7) 



where 



P(ri;t)=P + P (ri;*) + P (n; t) + P (n; t) , (8) 

g (n; f) = q k (n; f) + (n; f) + <z m/ (n; f) + g z (n; f) 

are the total stress tensor and vector of heat flow correspondingly. They have 
additive structure and contain several terms, each of them is stipulated by the 

influence from one of collision integrals [2,8]: P and q by Enskog collision 

integral (2), P and q mi by collision integral of KMFT (3), P and q by 

Landau collision integral (4), P and q are pure kinetic contributions only. 

P (ri,t) and q (ri,t) are new terms in the structure of (8) in comparison 
with results of [18]: 



P(r 1 ;t)= (9) 
Z 4 e 4 r , d f , r 12 r 12 r,o ( d d 



m 



r\ 2 g \ dvi dv 2 j J 



q l (ri;t)= (10) 

Z 4 e 4 r c2 _d_ r r 12 r 12 r 12 j d d) j dX p i 

2m J 1 1 dviJ 2 r\ 2 g \dvi dv 2 J J 

F l = <?2/i/i [8]. A short comment is needed for (7). First of all equation (7) is 
a balance equation for a kinetic part of total energy. To write the conservation 
law for total energy it is necessary to know also two-particle distribution func- 
tion f2{x\,x 2 ;t) next to one-particle one, because the potential part of total 
energy is expressed via f 2 . The Enskog-Landau kinetic equation in "pair col- 
lision" approximation has been obtained from the BBGKY hierarchy with a 
modified boundary condition in [2], where the expression for f 2 is also pointed 
out. An average value for the potential energy and its flow one should be cal- 
culated on the basis of this expression. Then, adding it to the balance equation 
(7), one can obtain the conservation law for total energy. 
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We shall construct a normal solution to the Enskog-Landau kinetic equation 
(1) using the boundary conditions method [12,13]. Following this method, let 
us bring into the right-hand side of equation (1) an infinity small source with 
e -> +0: 

Ie (xi, t) + I M f (xi; t) + I L (xi; t) - e (A (xi\ t) - / 1 (0) (xi, t)) , 

where A (0) ( Xl ;t) is some already known one-particle distribution function sat- 
isfying equations (5) - (7) for parameters of reduced description of our system. 
Then the solution can be found in the form /i (xi, t) = f[ ^ (x±; t) + Sf (x±; t) 
and search of the normal solution implies treatment of the correction Sf (x-y\t). 
Substituting fi(xi,t) into (11), one can obtain: 

{I + v 1 ± + s} Sf + JjLtf) = l MF (/(o)) + l MF {S f) + (12) 

Ie (A (0) , A (0)1 ) + h (A (0) , A (0) ) + Ie (f[°\ Sf) + I E (5f, f^) + 
h (A (0) , Sf) + I L (Sf, A (0) ) + I E (Sf, Sf) + I L (Sf, Sf) . 

Conventional signs used in the equation (12) are obvious [2,8,9]. Also the 
fact was taken into account about lMF(xi]t), collision integral (3), which is a 
functional of one-particle distribution function only. Terms with the subscript 
E are nonlocal, therefore in further calculations we should take their expansion 
with respect to the local one-particle distribution function and cut-off this 
expansion by terms with degrees higher than Sf. In the case when terms 
with subscripts MF and L also mean nonlocal functionals, one should apply 
mentioned above procedure to them too. Let us combine some terms in (12): 



Ie (Sf) = Ie Sf) + Ie (Sf, f[°^) linearized nonlocal Enskog colli- 
sion functional, 

II (Sf) = II (A^^ Sf) + I L (Sf, f[ ^) linearized Landau collision func- 
tional. 

Now let us designate L t (Sf) = 1^ (Sf) + Imf (Sf) + II (Sf) and introduce 
an operator S (t, t') with the following properties: 

(t, = L t (Sf) S (t, t>) , S (t, t') \ t , =t = 1. 
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Using these properties of operator S (t,t'), one can represent equation (12) 
in an integral form. Having correction Sf(x; t) in an integral form, it is easy 
to cross to itemizing procedure for finding it in corresponding approximation. 
For example, it can be organized in the following way: 



Sf (k+1) (*i;f) = / df e-^')S(M') l-§- t A 0) ~ ^^f ik) + (13) 

— oo * 

Ie (A (g) , A (0) ) + Imf (A (0) ) + h (a (0) , A (0) ) + i2W*)\ , 

where subscript t' at the bottom of right brase means that integrated expres- 
sion is a function of t'. An additional condition to find Sf(x;t) is the evident 
limit lim Sf(x; t) = 0. In order to construct the (fc + l)-th approximation it is 

necessary to use the fact that Sf\ k=0 = and the conservation laws (or equa- 
tions for reduced description parameters) in k-th approximation. To realize 
this procedure a zeroth approximation for the one-particle distribution func- 
tion fi°\xi,t) is needed. In the case of spherical charged particles, fi°\xi;t) 
can be chosen as the local-equilibrium Maxwell distribution function 

AO)/ +\ / +\( m / mci(ri,t)\ 



Let us find a correction to the distribution function t) using itemizing 

procedure (13). Calculating and obtaining of conservation laws (5), (6) and 
equation (7), we should take into account the following relations: 



92 (ri,r 2 ;*) = g 2 (r 1 ,r 2 ;n(t), P(t)) -> g 2 (r 12 ; n(t), /3(t)), 
F - F<°> = g 2 (r 12 ; n(t) , (3 (t)) fi° ] (*i5 t) A (0) (n, v 2 ; t) , 



(14) 



where g 2 (ri 2 ;n(t), (3(t)) is the binary quasiequilibrium correlation function, 
which depends on relative distance between particles. We obtain for stress 
tensor and heat flow vector: 



P k =7 P k , P k = nkT, (15) 

p hs =Y P hs , P hs = \ixn 2 a z kTg 2 (a\n, (3) , (16) 

oo 

p =lP mf i P"f = --n(nZe) 2 J—g2(r\n,P), (17) 
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o, 



,hs 



mf 



0. 



(18) 



«— ► ^ — Mi 

In these expressions / is the unit tensor, P and q are equal to zero because 
the integration between symmetrical limits goes over odd function. As far as 

calculated components P (ri;t), P (ri;t), P (ri;t) and P (r±;t) (15) - 
(18) are known, one can write total pressure in the zeroth approximation: 

/ 2 \ 2 °r dr 

P = nkT (l + -nna 3 g 2 (a\n,P)j - -ix{nZe) 2 J — g 2 {r\n,(3). 



Calculating expressions in brackets on the right hand side in (13), one can 
write total expression for correction 5f(x\\t) in first approximation: 



6fW ( Xl ;t) = - J dt' e- 6 **"* 7 ) S (t,t') 

— oo 

2 \ 
1 + -nn(r 3 g 2 {a\n,0)\ 



mc\ 5 
2kf ~ 2 



d 



(19) 



c 1 —\nT(r 1 ;t) + 
or i 



4 \ 777 

1 + Y5 nna3g2 ( a ' n '^)J ~tff 



1 



g 

9ri 



V(n;0 



Terms related to short-range interactions only contribute evidently into the 
correction in this approximation. Contrary to the kinetic theory of dilute gases 
particle sizes take part here [6,7,11], where particles are considered as point- 
like objects. Nevertheless, the influence of both long-range and short-range 
parts of interactions are also "hidden" in operator S(t, t') (through operator 
L t ). Formally, the expression (19) looks completely the same as the correction 
in [18]. But a difference lies in the structure of the operator S(t,t'). 

Having total expression for correction Sf(xi; t) in the first approximation (19) 
one can calculate conservation laws (5), (6) and equation (7) in the same ap- 
proximation. Therefore, it is necessary, first, to obtain relations for determining 

<->k 1 

quantities (7) in which the correction (19) can be engaged. For p (t*i; t) we 
obtain: 



p k 1 ( ri ; t) =7 P k - J dt' e - £{t '^ M k (t, t') 



(20) 



where S a p is a velocity shift tensor, 
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M k (t, *) = -] d Vl ddS (t, t') x 



[21) 
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A (0) (xi; t) (l + ^nn<7 S g 2 (a\n, /3)) ^ (c lCl - 7 
is a kernel of kinetic part of the transport equations. 

<->/is 1 <-ymf 1 <-+£ 1 

For calculating P (r*i; t), P (^i!*) an d -P (t*i; , we have to expand 
F hs , F m f and F l on inhomogeneity and deviation 5f (xi',t) and keep in the 
series initial terms only. The expansion for F m f , F l reads the same as for F hs 
with changing g 2 (a\n, (3) -> #2(1*12 K r 12 -> |r»i 2 |^Vi 2 , cr -> |r* 12 j - The 
calculations show: 



P (n; t) =P --n 2 a 4 g 2 (a\n, (3) VTimkT 
9 



6 <-> <-> 

+ (V • v / 

5 



(22) 



4 
15 



1 

7ina 3 g 2 (a\n, (3) J dt' e~ £ ^ M k (t, t') 



S 



(-tmf 1 <-> 

P ( ri ;t)=I P mf . 



(23) 



A mean field influence into the total stress tensor remains the same as in 



4 1 



zeroth approximation. Similar situation arises as to P (vi,t): 

P 1 (n;t) =P (n;t) = 0. (24) 



Total expression for stress tensor in the first approximation is a sum of (20), 
(22), (23) and (24): 



P (n;f) =/ - -n 2 a A g 2 (a\n, (3) Vnrnkf 

9 

t 

1 + — nna s g 2 (a\n, /3)) | dt' e^"^ M fc (t, t') 



6 



5+(V- V) / 



The calculations for heat flow vectors give: 



q k 1 ( ri; t ) = - J dt' e- £ ('-'') L k (t, t') 



T 



-VT 



—00 
2 
3' 



^ ul (r 1 ;t) = ~n 2 t T A g 2 (a\n, (3) 



lirk 3 T 



m 



vr(r i; t)- 



2 r 

-W^^n,/?) y dt' e^*') L fe (t, i') 



(25) 
(26) 



8 



g ,1 (r 1 ;*) = q , (r 1 ;*)=0. (27) 

Here 



(28) 



is another kernel of kinetic part of transport equations. Total expression for 
heat flux vector is a sum of (25) - (27): 



q(ri,t) = --n 2 a 4 g 2 (a\n, f3) 



I nk 3 T 



m 



VT(n;t)- 



1 + -irna 3 g 2 (a\n, /3)) J dt' e'^'^ L k (t, t') 



T 



-VT 



Now we can consider one of the limiting cases, namely, the stationary process, 
when the operator L t does not depend on time, i.e. S(t,t') = exp {L t (t — t')}. 

Some terms in expressions for p (t*i; t) and q{r\] t) can acquire simpler form. 
We can compare them with those from the Enskog-Landau kinetic equation 
for one-component system of charged hard spheres with using the Chapman- 
Enskog method in the case, when in a long-range part of the collision inte- 
gral we put g 2 (a\n, f3) — > 1. It should be noted that bulk viscosity has the 
same structure as in the Chapman-Enskog method [8]. But other transport 
coefficients exhibit some distinctions. The structure for shear viscosity rj and 
thermal conductivity A is: 



4 .3, 



,(<T|n,/3)J 
5 ~ ' {l ( E\Sf) + I L (5f)} 



3 < i + Tzima g 2 { 
v = - as + 2 nkT 1 f , x f (29) 



f 2 1 2 

U \ 1 + -7ma 3 g 2 (a\n, (3) \ 

2m m {4°)(5/) + / L (5/)} 

Then the problem lies in calculating collision integrals 1^(5 f) and I^Sf), 
this means that we should calculate collision integrals (2) (in the zeroth ap- 
proximation on inhomogeneity) and (4) together in the first approximation 
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on deviation 5f, where Sf is substituted from (19). The matter of some diffi- 
culty is that correction (19) in its turn is expressed also via collision integrals 
Ifi\5f), I^Sf), which are in the operator S(t,t'). So the first acceptable ap- 
proximation should be that, when correction 5f (19) is expressed via Ig\Sf), 
I L (Sf) calculated with 5f, where 5f = 5f at S(t,f) = 1. For if^f) we 
obtain the results [8], for ii,(<5/) in (29), (30) we can obtain the following: 



h{of) = — YJT~ / dri2 dv2 92 ri + r 12 ; t — g ^ — S 

at?i J rj 2 g [cwi cw 2 J 

x {/i (xi; 5f (n + ri 2 , v 2 ; i) + 5/ (x i; *) /i (n + r 12 , u 2 ; t) }, (31) 
where 5/ (x; t) is evaluated from (19) with 

S(t,t') = exp {L t (t - 0} = exp {[/^V/') + - f)} 



at Sf (x;t) = 5f(x;t) <^ ^ _ ^. This stage of calculations needs an ex- 
plicit form for the binary quasiequilibrium correlation function g 2 both on the 
contact and in r-space. 

The results of this paper (19), (29) and (30) will coincide completely with 
those from [8] when in a long-range part of the collision integral (31) one 
puts ^2(^1,^1 + r i2) = 1 and represents it in Boltzmann-like form. But the 
used boundary conditions method has turned out more convenient than the 
Chapman-Enskog one [2,8]. As was discussed in details in [12] at construct- 
ing the normal solution for a kinetic equation using the boundary conditions 
method, time derivatives d/dt of hydrodynamic parameters of reduced descrip- 
tion do not set to be small. Therefore, the normal solution to this equation 
could be used for hydrodynamic description of fast processes. 
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